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Abstract 

In this paper we explicitly work out the effects that a spherically 
symmetric distribution of dark matter with constant density would in- 
duce on the Keplerian orbital elements of the Solar System planets and 
compare them with the latest results in planetary orbit determination 
from the EPM2004 ephemerides. It turns out that the longitudes of 
perihelia w and the mean longitudes A are affected by secular preces- 
sions. The resulting upper bounds on dark matter density, obtained 
from the EPM2004 formal errors in the determined mean longitude 
shifts over 90 years, lie in the range 10~ 19 — 10~ 20 g cm -3 with a peak 
of 1CP 22 g cm -3 for Mars. Suitable combinations of the planetary 
mean longitudes and perihelia, which cancel out the aliasing impact of 
some of the unmodelled or mismodelled forces of the dynamical mod- 
els of EPM2004, yield a global upper bound of 7 x 10~ 20 g cm~ 3 and 
4 x 10~ 19 g cm -3 , respectively. 
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1 Introduction 



In this paper we exploit the latest accurate data on planetary orbit deter- 
mination in order to constraint the amount of dark matter present in 
the Solar System. To avoid possible misunderstandings, here we are not 
intending the amount of dark matter at the Solar System's position in the 
Galaxy; we are, instead, dealing with possible excess of dark matter above 
the galactic background density 1 . Such excess is to be intended as coming 

1 Such a galactic background density would not cause the perturbing effects on planetary 
motions discussed in this paper; it is believed to cause an overall acceleration toward the 
center of the Galaxy. 
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from dark matter gravitationally bound to the Solar System and partly from 
dark matter in the galactic halo focused by the solar gravitational potential. 

Let us consider a distribution with spherically symmetric density pdm(r)- 
The corresponding modification to the Newtonian gravitational potential 
can be cast in the form 

U dm = -^-, (1) 
r 

where 

H(r) = 4vr / p d m(r')r' 2 dr' (2) 
Jo 

is the total mass of dark matter enclosed in a spherical region of radius r. 
In regard to the distribution of the dark matter, it is believed it comes from 
diffusion into the Solar System (or possibly it was already present when our 
planetary system was formed). Diffusion was studied, e.g. in jHO], but it 
turns out that in none of these papers the density distribution in the Solar 
System has been derived in a form suitable for the present paper. Thus, we 
will make the simplest choice, i.e. we will consider pd m as constant. With 
this assumption p = (A/3)npd m r s and eq. (^Q) becomes 

Udm = ~^Gp dm r 2 . (3) 

From it an entirely radial acceleration 

g 

-4dm = ^Gpdmr (4) 

occurs. Its effect on planetary motions, which is, of course, much smaller 
than usual Newtonian gravity, can straightforwardly be calculated within 
the usual perturbative schemes. 



2 The impact of a spherically symmetric distribu- 
tion of dark matter on the planetary orbits 

The Gauss equations for the variations of the semimajor axis a, the ec- 
centricity e, the inclination i, the longitude of the ascending node O, the 
argument of pericentre u> and the mean anomaly M of a test particle in the 
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gravitational field of a body M are 
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in which n = 2ir/P is the mean motion 2 , P is the test particle's orbital 
period, / is the true anomaly counted from the pericentre, p = a(l — e 2 ) 
is the semilactus rectum of the Keplerian ellipse, A r , At, A n are the ra- 
dial, transverse (in-plane components) and the normal (out-of-plane com- 
ponent) projections of the perturbing acceleration A, respectively, on the 
frame {f,t,h} comoving with the particle. 

If the perturbing acceleration is entirely radial, as it is the case for eq. 
(UJ), the Gauss equations reduce to 
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As a result, the inclination and the node are not perturbed by such a 
perturbing acceleration. 



2 For an unperturbed Keplerian ellipse it amounts to n = GM/o 
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By evaluating eq. @ on the unperturbed Keplerian ellipse 



^ - e2 ) (17) 
1 + ecos/' [ ' 



inserting it into (|lll) - (|16|) and averaging them with 

dt 1(1- e 2 ) 3 / 2 



P 2vr (1 + ecos/)< 



df (18) 



one gets that only the argument of pericentre and the mean anomaly are 
affected by extra secular precessions 

duo \ AirGpdr, 



Vl-e 2 , (19) 



dt I n 
iM\_^ {20) 



As a consequence, also the mean longitude A = Q + uj + M, which is used for 
orbits with small inclinations and eccentricities as those of the Solar System 
planets, is affected by an extra secular rate 



3 Confrontation with the latest data 
3.1 The longitudes of perihelia 

The obtained expressions are useful for comparison with the latest data on 
planetary orbits from the EPM2004 ephemerides . A similar analysis was 
conducted in [I] by exploiting the determined extra-advances Azu of the 
longitudes of perihelia w of the inner planets By using the uncertainties 
5(Aw) released in 0, the authors of [3] obtained an upper bound on p^m 
of 10 -17 g cm -3 from Mercury and of 10~ 19 g cm~ 3 from the Earth and 
Mars. Note that eq. ()19|) yields for the perihelion advance, in units of 2tt, 
after one full orbital revolution the same expression (5) of (3], apart from 
the eccentricity factor. Moreover, the constraints on p^ m obtainable for each 
planet separately with eq. (fT9|) and the results of [2] fully agree with those 
of [3]. However, the following concern about the reliability of the estimates 
by the authors of [3J is in order. The extra-perihelion advances determined 
in [Sj are affected, in general, by all the Newtonian and non-Newtonian 
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features of motion which have not been accounted for in the dynamical force 
models of EPM2004. Among them there are certainly the totally umodelled 
general relativistic gravitomagnetic field, which induces the Lense-Thirring 
planetary precessions, and the solar quadrupole mass moment J2 which, 
instead, was included in EPM2004, but it is currently affected by a ~ 10% 
uncertainty. For a discussion of such issues in a different context see [JJ]. 
Each planetary perihelion is, thus, affected by such residual effects, so that 
it is not possible to entirely attribute the determined extra-advances to the 
action of the dark matter distribution, especially for Mercury. A better 
way to use Am is to suitably combine them in order to make the estimate 
of pdm independent, by construction, of the Lense-Thirring and J2 effects. 
According to the approach followed in [SJ [7| , it is possible to construct the 
following combination 

Ati7 Mer + ciAro Ear + c 2 ACT Mar 
Aim - 1.7 x i io cm 3 y,.-! g -l ' W 

with ci = —80.7 and C2 = 217.6. The dimensionless coefficients c\ and 
C2, which are built up with a,e and i of the planets adopted, cancel out, 
by construction, the impact of the gravitomagnetic field and of the solar 
quadrupolar mass moment on the combination eq. (|22[). This can straight- 
forwardly be checked by combining the perihelion precessions due to J2 and 
to the Lense-Thirring force with the coefficients of eq. (|22jl: the result is 
zero. The upper bound on p^m obtainable from eq. (|22|) can conservatively 
be obtained from 3 

A . 8(Am Me <) + \ Cl \5(Aw E ^) + \c 2 \5(Am M ™) 
Pdm ~ 1.7 x 10 10 cm 3 yr- 1 g- 1 ' { 6) 

it amounts to 4 x 10~ 19 g cm -3 . 



3.2 The mean longitudes 

Here we compare the shifts in A obtainable from eq. (|21|) with the formal 
errors in fitting the mean longitudes at epoch of all the nine major planets 
by using a data record spanning 90 years from Table 4 of [Q . Our results are 
in Tabled As can be noted, the tightest constraint comes from Mars, with 
5 x 10" 22 g cm~ 3 . The present upper bounds are better than those obtained 

3 Although rather small, a certain amount of correlation among the determined extra- 
advances of perihelia is present, with a maximum of 20% between the Earth and Mars (E.V. 
Pitjeva, private communication, 2006). If a root-sum-square calculation is performed, an 
upper bound pdm < 3 x 10~ 19 g cm" 3 is obtained. 
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Table 1: Upper bounds, in g cm -3 , on p^ m from the comparison of the 
planetary mean longitude shifts over 90 years from eq. (|21|) and the formal 
errors in the planetary mean longitudes released in PQ. 



Mercury 


4 x 10- 


-19 


Venus 


9 x 10- 


-20 


Mars 


5 x 10- 


-22 


Jupiter 


2 x 10- 


-20 


Saturn 


3 x 10- 


-20 


Uranus 


3 x 10" 


-20 


Neptune 


6 x 10- 


-20 


Pluto 


9 x 10" 


-20 



Table 2: Secular precessions of longitudes of perihelia and mean longitudes, 
in arcseconds per century, induced on the outer planets by a dark matter 
density p^ m = 10~ 16 g cm -3 equal to the upper bound found in jSj. 



Planet 


w 


A 


Jupiter 


3.25 


-4.349 


Saturn 


8.06 


-10.803 


Uranus 


23.0 


-30.807 


Neptune 


45 


-60.27 


Pluto 


65.8 


-97.03 



in [I] from the perihelia of each planet, even if we conservatively assume that 
the realistic errors in A may be one order of magnitude larger than those 
released in pp. In regard to the constraint from Uranus 3 x 10 -20 g cm~ 3 , 
the previously obtained upper bound on p^ m from such a planet amounted 
to 10 -16 g cm -3 jSj. From Table 121 it can be noted that with a such high 
value of pdm the periehelia and the mean longitudes of the outer planets 
of Solar System would undergo very huge secular precessions up to tens of 
arcseconds per century. Such anomalous effects would have certainly been 
detected, even in the case of the outer planets whose orbital determination 
relies almost entirely upon relatively inaccurate optical data [Q. 

However, caution is advised in considering the results of Table^for each 
planet separately. It is so not only for the same reasons already exposed in 
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Section |3~T1 about the unmodelled or mismodelled forces, but also because of 
the systematic errors 5n = (3/2) ^jGMja^ba in the Keplerian mean motions 
induced by the uncertainties in the semimajor axes 4 . Following the approach 
of Section 13.11 we can write down the combination 

_ AA Mcr + fciAA Vcn + fc 2 AA Mar 
Pdm ~ -2.4 x 10 10 cm 3 yr" 1 g" 1 ' { ' 

with k\ = —15.4 and k,2 = 78.1. Also such a combination is independent of 
the solar Lense-Thirring and J2 effects. The formal errors in A of Table 4 
of pQ, combined according to eq. (|24jl. yields a conservative upper bound 
Spdm < 7 x 10~ 21 g cm -3 . The errors in n, linearly combined according to 
eq. (|24j). yield a systematic bias 5pd m < 7.7 x 1CP 20 g cm -3 , so that a total 
upper bound 5pd m < 7.8 x 1CP 20 g cm -3 can be obtained with the squared 
sum of these two kinds of errors. As can be noted, the use of eq. (01)) yields 
better results that eq. (|22j). 



4 Conclusions 

In this paper we have worked out the effects that a local excess of dark mat- 
ter in our Solar System over the galactic background would induce on the 
orbits of the planets. It turns out that, by modelling the dark matter dis- 
tribution with a constant density, the longitude of perihelion and the mean 
longitude of a planet are affected by secular precessions. The comparison 
with the latest data show that the upper bounds obtainable from the mean 
longitudes and the perihelia are of the order of 10 -20 g cnr 3 and 10~ 19 g 
cm~ 3 , respectively. 
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